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Introduction 

The  Xa  method1  and  its  offspring  are  becoming  more  and  more  im¬ 
portant  in  quantum  chemistry  for  chemical  systems  that  are  too  large  for 
accurate  configuration  interaction  (Cl)  calculations.  All  self-consistent-field 
(SCF)  methods  of  quantum  chemistry  method  begin  with  one  or  more  sys¬ 
tems  of  one-electron  equations  of  the  form, 

=  [-iV2  +  l/]*.(r).  (1) 

Invariably  V  is  separated  into  its  one-electron  and  two-electron  parts,  the 
electron-nuclear  interaction  potential,  Ven,  and  the  electron-electron  inter¬ 
action  potential,  Vee.  respectively.  Only  Vtf.  distinguishes  the  various  SCF 
methods  of  quantum  chemistry.  In  what  are  now  called  local  density  func¬ 
tional  (LDF)  methods  the  electron-electron  potential  is  written, 

vtt  =  V«(ri)  =  [  f,{y-}-V-  +  V,c(p(n)),  (2) 

where  the  total  electron  density,  p.  can  be.  divided  into  spin-up  and  spin- 
down  components,  and  the  expression  for  the  spin-up  component, 

o(1-)  =  », ,  <?;,(r)<pi,  (r),  (3) 

where  the  rq  are  the  (possibly  fractional)  number  of  electrons  occupying 
each  orbital  of  Eq.  1,  contains  an  implicit  summation  over  orbital  index 
i.  These  methods  are  called  LDF's  because  the  exchange  and  correlation 
potential,  VIC(p{i'\ )),  at  any  point  in  space  is  constrained  to  depend  only 
on  the  density  at  the  same  point.  If  we  constrain  VIC  to  be  local  and 
require  the  system  of  equations  to  satisfy  the  virial  theorem  at  all  extrema 
on  all  I3orn-Oppenheimer  potential  energy  surfaces,  we  are  led  to  the  Xo 
approximation  for  spin-up  orbitals. 

V'AtMf.))  =  3o(^p,(r))‘/3.  (4) 
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This  result  follows  from  simply  considering  the  homogeneous  electron  gas 
at  all  densities  and  requires  that  LDF  expressions  not  include  gradients  of 
the  density.  The  value  of  a  for  the  homogeneous  electron  gas  based  on 
variation  is  two-thirds2  and  based  on  averaging  is  one.1  For  heteronuclear 
molecules  the  best  choice  is  perhaps  0.7. 

On  the  other  hand,  the  Kohn-Sham2  mapping  of  the  the  ground  state  of 
the  real  system  onto  the  ground  state  of  a  fictitious  noninteracting  set  of 
electrons,  which  has  the  same  density  as  the  real  system,  provides  a  different 
interpretation  of  these  equations.  If  both  the  real  and  noninteracting  sets 
of  electrons  are  nondegenerate,  the  mapping  is  one-to-one  and  the  wave- 
function  of  the  noninteracting  set  of  electrons  is  a  single  determinant  of 
one-electron  orbitals  likely  satisfying  Eq.  1  for  some  V .  (The  counterexam¬ 
ples  of  non-'Z-iepresentability  use  degeneracy  ol  one-electron  orbitals. ■5~s) 
If  we  map  the  energy  of  a  single-determinant  of  plane-wave  orbitals  onto 
the  essentially  exact  energy6  of  the  homogeneous  electron  gas  as  a  function 
of  background  positive  charge  density  in  the  completely  spin-paired  and 
completely  spin-polarized  extremes,  we  are  led  to  a  different  set  of  LDF’s 
that  typically  overbind  molecules  compared  to  Xa  and  experiment.7  The 
best  of  this  other  set  is  perhaps  the  Perdew-Zunger  (PZ)  parameterization,8 
which,  although  it  is  too  complicated  to  write  out  here,  gives,  apart  from 
slightly  rescaling  the  total  energies,  (  lie  same  description  of  most  nonmag¬ 
netic  molecules  as  does  Xo.  The  reason  for  almost  no  difference  is  that 
the  one-third  power  and  functions  similar  to  it  tremendously  deamplify 
variations  in  the  nonnegative  density  itself 

The  computational  attraction  of  Eq.  1  and  2  is  that  V'(r)  is  a  local  po¬ 
tential,  i.e.,  it  is  the  same  for  each  orbital. While  all  the  orbital  generating 
equations  of  quantum  chemistry  can  be  rendered  in  the  form  of  Eq.  I  (if  all 
else  fails  by  using  projection)  only  with  a  local  potential  is  l'(r)  a  single¬ 
valued  (orbital  independent)  function  of  position.  Thus  obtaining  V'(r) 
is  central  to  all  local  potential  methods  and  that  problem  is  straightfor¬ 
ward.  Eqs.  1-3  can  be  solved  for  any  real  vector  of  occupation  numbers, 
n.  All  LDF  approximations  to  l Ic  yield  local  potentials,  but  not  all  local 
potentials  result  from  LDF  approximations,  as  is  obvious  from  the  direct 
Coulomb  term  in  Eq.  2. 

For  any  density  functional  expression  for  lrf  the  coupled  one-electron 
equations  of  motion.  Eq.  1.  can  be  integrated  to  yield  the  total  energy 
expression, 


i:  =  (T)  +  (U), 


(5) 


where  the  total  kinetic  energy  is  given  by  an  expression, 
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(T)  =  -i«i,0;,(r)V2^<(r),  (6) 

containing  an  implicit  summation  over  the  two  spin  directions  through  the 
index  s,  and 


{U)  =  J  p3(v)Us(p(v))dr,  (7) 

wliere  the  one-electron  potential  and  the  density  functional  potential  energy 
operator,  U ,  are  related  by  an  independent  variation  of  (U)  with  respect  to 
each  spin  density  at  each  point, 

J  V^(p)6P](v)dr  ~  6(U)  (8) 

For  any  part  of  V  that  is  approximated  as  being  an  LDF,  Eq.  8  is  precisely 
a  differentiation, 

=  —  (l')^;(!)('’(r))l  Pi 

That  Eq.  8  does  not  hold  for  a  nonLDF  contribution  to  V  is  obvious  if  cne 
considers  multiplying  Eq.  2  by  /j(iq).  In  that  case  the  p( iq)  and  p( r2)  in 
the  direct  Coulomb  term  are  asymctrical;  r2  is  an  integration  variable  while 
ri  is  not. 

As  is  the  case  with  Hartree-Fock- based  quantum  chemical  methods,  the 
analytic  basis  set  approach  is  the  overwhelming  favorite  in  nonmuffin-tin 
local  potential  quantum  chemical  methods,  except  for  lower  dimensional 
problems  such  as  atoms10  and  linear  molecules.11-13  The  all-purpose  an¬ 
alytic  basis  set  is  Gaussians'  14  in  both  quantum  chemical  methods.  In 
contrast  to  Hartree-Fock  based  methods,  computational  efficiency  dictates 
and  complicated  expressions  such  as  Eq.  -1  inspire  fitting  V  in  analytic-basis- 
set  local  potential  methods.15-1''  If  the  potential  is  fit,  which  is  indicated 
by  placing  a  bar  over  the  fitted  quantity,  then  the  potential  energy  should 
be  evaluated,17 


(U)  «  J  p(r)V(P( r))dv  +  j  p(v)[U(p(v))  -  r(Mi-))]c/r,  (10) 

because  then  the  nonfitted  density  only  occurs  multiplied  by  the  local  poten¬ 
tial.  Variation  with  respect  to  t lie  orbitals  (occurring  only  in  the  nonfitted 
density)  gives  the  local  potential  term  of  the  one-electron  equations,  and 
the  total  energy  is  insensitive  in  first  order  to  changes  of  the  one-electron 
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orbitals,  provided  the  fitting  procedure  follows  variationally  from  this  equa¬ 
tion.  Eq.  10  yields  an  energy-stationary  way  to  fit  the  charge  density,7-18  to 
fit  the  direct  electron-electron  Coulomb  potential  itself,19  and  to  fit  the  Xa 
exchange  potential.7,20  For  the  more  complicated  Vzc  such  as  occur  in  the 
PZ  functional,  Eq.  10  yields  variational  weights  for  any  numerical  fitting 
scheme.21 

Independent  of  how  V  is  determined,  it  is  natural  to  ask  what  symmetry 
it  has  and  to  what  extent  symmetry  can  help  in  local  potential  calcula¬ 
tions.  The  system  of  equations  is  nonlinear,  and  therefore  its  solutions  can 
and  do  yield  densities  that  break  the  symmetry  of  the  collection  of  nuclear 
charges.22  There  are  three  approaches  to  this  important  problem.  First, 
one  can  accept  broken-symmetry  solutions.23-26  Second,  one  can  restrict 
one’s  attention  to  density  functional  methods  that  yield  V  having  the  sym¬ 
metry  of  the  nuclear  charges.3-5  Third,  one  can  explicitly  symmetrize  V 
by  symmetrizing  a  broken-symmetry  wavefunction  before  computing  the 
kinetic  energy  and  density  during  the  SCF  process.27 

Empirically,  the  experimental  ionization  potential  of  a  core-hole  from  a 
set  of  symmetry-equivalent  atoms  is  better  reproduced  in  local  potential 
methods  when  symmetry  is  broken  and  the  hole  self-consistently  localizes 
on  one  of  that  set  of  atoms.23,24  Unfortunately,  in  contrast  to  ab  initio 
methods,  the  variational  principle  does  not  also  suggest  this  solution.  In 
current  local  potential  methods,  a  localized  core  hole  results  in  a  higher 
total  energy  than  a  delocalized-hole  total  energy.  When  applied  to  spin, 
symmetry-restriction  leads  to  a  poor  description  of  Cr2, 25,26  which  has  a 
singlet  electronic  ground  state  but  must  dissociate  into  two  septet  atoms. 
All  current  LDF  methods,  however,  are  not  invariant  under  rotations  in 
spin  space.22,26  Therefore  they  can  never  yield  spin-densities  that  trans¬ 
form  properly  for  any  magnetic  system  that  has  any  minority-spin  elec¬ 
trons,  without  further  assumptions  (such  as  using  the  same  orbitals  for  the 
two  spins).  Nevertheless,  significant  progress  is  being  made  in  interpreting 
solutions  in  which  the  density  is  invariant  but  the  spin  densities  break  the 
molecular  symmetry29,30  (spin  density  waves).  Apart  from  the  special  case 
of  core  holes,  there  has  been  no  general  attempt  to  interpret  broken-spatial- 
symmetry  solutions  (charge  density  waves)  despite  the  fact  that  they  they 
often  arise  asymptotically  in  molecular  dissociation.11,31  In  cases  where  the 
real  system  does  not  exhibit  cliarge-density-wave  behavior,  these  broken- 
spatial-symmetry  errors  typically  make  binding  energies  uncertain  at  the 
tenth  of  an  electron  volt  level. 

Symmetrizing  V  removes  this  uncertainty.  (For  atoms  this  means  mak¬ 
ing  the  central  field  approximation.)  Furthermore,  a  symmetrized  V 
can  be  readily  taken  advantage  of  in  linear-combination-of-atomic-orbitals 
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(LCAO)  local  potential  computer  codes.22  Symmetpzation  without  further 
approximation  allows  routine  LDF  calculations  on  symmetric  systems  con¬ 
taining  tens  of  transition  metal  atoms,32-3,4  or  larger  mixed  system  such  as 
C6o  and  yttrium  diphtalocyanine  (YCe-il^Nis).35'36 

In  what  follows  a  density  functional  theory2''26  of  the  broken-spatial- 
symmetry  problem  for  singlet  electronic  states  is  reviewed.  It  relies  on 
the  fact  the  density  is  to  be  symmetrized,  which,  in  turn,  symmetrizes  V 
provided  it  is  a  function  of  the  density.  This  symmetrization  validates  the 
fractional  occupation  number  (FON)  description  of  the  ground  state.37  The 
theory  contains  the  mathematics  of  ab  initio  configuration  interaction.  A 
good  introduction  to  any  discussion  of  correlation  and  LDF  theory  is  the 
work  of  Cook  and  Karplus.36 

FON  Singlet-State  Configuration  Interaction 

The  simplest  densities  to  consider  correspond  to  closed-shell  molecules, 
in  which  the  density  clearly  has  the  symmetry  of  the  molecule.  For  such 
a  molecule  consider  the  case  where  the  highest  occupied  molecular  or¬ 
bital  (HOMO)  is  antibonding  and  the  lowest-unoccupied  molecular  orbital 
(LUMO)  is  bonding  and  compress  the  appropriate  internuclear  separations 
until  these  orbitals  become  nearly  degenerate  (or  consider  the  opposite 
case  and  expand  the  bonds).  For  concreteness,  consider  a  4E+  state  of  a 
homonuclear  diatomic  molecule  and  a  broken-symmetry  HOMO  of  mixed 
d g  and  <ju  character, 


4>  =  a\o,  +  Ho.  (U) 

The  g  and  u  components  must  be  orthogonal  from  symmetry  considerations 
alone.  (If  they  are  also  normalized,  then  the  sum  of  the  magnitudes  squared 
of  a  and  6  is  unity.)  Since  the  kinetic  energy  operator  is  invariant  under 
any  spatial  symmetry  operation,  o’s  contribution  to  the  kinetic  energy, 

T*  -  -a-a(xvf|iV2|.V>,)  -  b'b  {Xa.  |iV2|.X,.),  (12) 

is  diagonal  in  its  components;  thus  the  kinetic  energy  is  the  same  as  that 
of  asymmetry-restricted  FON  calculation  using  na  -  a'o  and  nff,  =  b'b. 
In  contrast  to  the  situation  for  the  overlap  and  kinetic  energy,  there  is  a 
difference  between  the  broken-symmetry  and  symmetry-restricted  densities. 
The  difference  is  the  unsymmetrica!  expression, 

a‘b  uo„uo.  +  b’a  (13) 

that  has  larger  magnitude  on  either  the  left-hand  atom  or  on  the  right 
hand  atom  depending  on  the  magnitudes  and  phases  of  a  and  b.  If,  how¬ 
ever,  the  density  is  symmetrized  by,  in  this  case,  averaging  the  original 


6 


broken-symmetry  orbital  density  with  the  same  density  after  inverting  the 
coordinate  system,  then  botli  orbital  densities  are  the  same, 


a'“  u'°,uo,  +  b'b  (14) 

Thus  symmetrized-broken-symmetry  and  the  FON  density-functional  cal¬ 
culations  are  identical.  For  the  above  example,  both  give  a  charge  density 
consistent  with  a  1 E*  electronic  state. 

FON  SCF  calculations  are  always  slower  to  converge  than  fixed-occupa- 
tion-number  LDF  calculations.  This  is  precisely  because  there  is  at  least 
one  more  degree  of  freedom  involved.  Another  way  to  analyze  this  slow¬ 
ness  is  to  note  that  for  fixed  occupation-number  calculations,  the  relevant 
gauge  of  speed  of  convergence  is  the  gap  between  the  HOMO  and  LUMO 
energies,39  whereas  FON  convergence  is  also  strongly  influenced  by  the 
density  of  states  surrounding  the  HOMO  and  LUMO  energies.17  (For  the 
FON  case,  the  HOMO  must  be  defined  as  the  highest  level  that  is  greater 
than  half  occupied  and  the  LUMO  is  the  lowest  level  that  is  less  than  half 
occupied.)  Therefore  it  is  best  to  avoid  direct  FON  calculations  if  one  can. 

To  an  excellent  approximation  the  total  LDF  energy  is  a  quadratic  func¬ 
tion  of  the  occupation  numbers,37  because  all  of  its  components  except  Vxc 
are,  and  Vrc  is  a  very  slowly  varying  function  of  the  density.  This  quadratic 
property  is  also  suggested  by  the  fact  that  the  exact  total  energy  depends 
only  on  the  first  and  second  order  density  matrices.  If  the  LDF  energy  is  a 
quadratic  function  of  the  occupation  numbers,  then  the  LDF  one-electron 
eigenvalues,  which  are  the  derivatives  of  the  total  energy  with  respect  to 
the  corresponding  occupation  number,40 


dE 

On,, 


(15) 


must  be  linear  functions  of  the  occupation  numbers.  In  particular,  this 
approximation  means  that  for  the  a9  and  cru  orbitals  under  consideration, 


At  =  <u  -  c9  (16) 

varies  linearly  with  fractional  number  of  electrons,  n,  transferred  from  the 
a g  orbital  to  the  cru  orbital.  The  value  of  this  eigenvalue  difference  for  all 
numbers  of  electrons  transferred  then  follows  from  the  lagrange  interpolat¬ 
ing  formula, 


Ac(n)  =  At(0) 


(n  ~  2) 
(0-2) 


-  Ac(2) 


(n-0) 

(2-0)' 


(17) 


which  reduces, 


(18) 


M»)  =  At(0)  -  [(Af(2)  —  Af(0)] 
There  is  an  FON  solution  A e(ny)  =  0,  if  and  only  if 


Af(0)  Af(2)  <  0.  (19) 

In  that  case,  transferring  nj  electrons. 


2Ac(0) 

W/  ~  Ac(0)  -  Ae(2)' 


(20) 


gives  the  FON  ground  state. 

Integrating  to  get  the  total  energy  at  the  FON  solution  from  the  n  =  0 
limit, 


f"' 

A  E  =  /  A  e(n)dn 

Jo 

gives  one  formula  for  the  FON  total  energy, 


Ej  o  =  £(0)  + 


Af(0)2 

Ac(0)  -  Ac(2) ’ 


and  integrating  from  the  other  limit  gives  another  formula, 


(21) 


(22) 


E, ,  =  £(  2)  + 


Ae(2)2 

Ac(0)  -  Ac(2) 


(23) 


Lagrange  interpolating  between  these  two  expressions  for  the  FON  total 
energy  gives  an  FON  energy,26 


— £(0)Ac(2)  +  E(2)A<(0)  -  Ae(0)Ac(2) 
Af(0)  -  Af(2) 


(24) 


that  connects  continuously  to  the  pure  state  solutions  outside  the  FON 
range,  geometries  for  which  Eq.  19  is  not  satisfied. 

The  FON  energy  lowering  can  be  viewed  as  resulting  from  configuration 
interaction  between  the  ;»  =  0  and  the  it  =  2  states,  in  which  the  upper 
FON  energy, 


E'(0)Af(0)  -  g(2)Ac(2)  +  Ac(0)Af(2) 
A<(0)  -  Ac(2) 


(25) 


lies  as  much  higher  in  energy  above  the  average  of  E{ 0)  and  £7(2)  as  Ei 
lies  below  that  average.  This  approach  can  be  extended  to  the  general 
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n  x  n  configuration  interaction  problem.  In  such  an  extension  the  pure 
state  energies  become  diagonal  entries  in  an  n  x  n  eigenvalue  problem.  In 
this  eigenvalue  problem  one  adds  an  ofif-diagonal  matrix  between  any  two 
pure  states, 


if  and  only  if  the  two  pure  states  differ  precisely  in  that  two  electrons  are 
transferred  from  one  orbital  to  another  and  that  the  eigenvalue  difference 
between  these  two  orbitals  change  according  to  Eq.  19  in  going  from  the 
one  state  to  the  other.  This  postulate  gives  as  its  ground  state  the  FON 
solution  when  restricted  to  only  allowing  a  single  pair,  but  any  single  pair, 
of  states  to  exchange  electrons. 

AI4 

The  Xa  ground  state  of  AI4  is  a  D-^/, -symmetric  rhombus  with  triplet 
electronic  configuration  having  a  spin-up  electron  in  each  of  the  9aIy  and 
16i„  orbitals  outside  a  50-electron  8^26^36^  closed-shell 

configuration.  The  big  orbital  is  a  x  orbital  with  nodal  plane  in  the  plane 
of  the  molecule.  This  agrees  with  an  <16  initio  description. 'u  The  triplet- 
state  atomization  energy  (the  energy  required  to  dissociate  the  molecule 
into  four  separated  R(3)-synmietrizcd  spin-polarized  atoms)  is  6.52  eV,  its 
bond  distance  is  4.99  bolir,  and  it  has  a  bond  angle  of  88°. 

Much  more  interesting  than  the  ground  state  potential  energy  surface  are 
the  singlet  potential  energy  surfaces.  Fig.  la  gives  the  symmetry-restricted 
low-lying  surfaces  in  square-planar  D-w,  symmetry.  These  breathing-mode 
surfaces  are  derived  from  the  51-clcctron  lu]glalg'2a2u4b]g‘2blslbnulegbc* 
configuration  by  anhilating  two  electrons  from  the  valence  a\g,  02U,  6js, 
and  bn g  orbitals  respectively.  There  arc  five  curve  crossings  in  the  figure. 
This  case  is  ideal  for  using  Eqs.  24-26  in  a  4  x  4  eigenvalue  problem.  The 
result  is  the  solid  lines  of  Fig.  1  b,  where  all  crossings,  of  the  original,  dashed 
lines,  are  avoided. 

The  most  interesting  avoided  crossing  in  Fig.  lb  is  the  one,  at  smallest 
Al-Al  bond  distance,  between  the  n°v  and  states.  It  is  the  weakest.  This 
is  because  the  aig  orbital,  being  nodeless,  is  largely  centered  on  the  atomic 
centers,  like  the  b\g  orbital  which,  being  t/x3_y:-like,  has  nodes  at  the  bond 
centers  for  atoms  chosen  as  they  were  to  lie  in  the  ±z  and  ±y  directions. 
Because  the  densities  ar^  similar,  transferring  two  electrons  between  these 
two  orbitals  has  the  smallest  effect  on  their  eigenvalue  difference,  and  the 
configurations  do  not  repel  each  other  very  much. 

This  is  not  the  whole  story  for  singlet  AI4,  by  any  means.  Fig.  2a  consid¬ 
ers  D-j/i-syimnetric  bending  motion  at  the  minimum  Al-Al  bond  distance 
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of  4.9  boiir  in  Fig.  1,  where  the  ground  state  is  largely  a\g.  Under  this  mo¬ 
tion,  the  a\g  state  interacts  with  two  ;>like  states  that  have  nodes  along  the 
shortest  Al-Al  second  neighbor  bond  when  they  cross  the  a\g  state.  These 
two  p-like  states  are  degenerate  in  symmetry  and  thus  are  different, 
higher  lying  in  energy,  from  the  states  considered  in  Fig.  1.  Fig.  2b  gives 
the  result  of  using  Eqs.  24-26  on  the  pure  states  of  Fig.  2a. 

Conclusions 

The  work  of  def.ning  local  potential,  Xa-like,  methods  for  use  in  quantum 
chemisty  is  not  done.  Only  pure  singlet  states  were  addressed  in  this  work, 
and  only  methods  that  resulted  in  no  net  spin-polarization  anywhere  in 
space  is  space  were  used.  The  problem  of  coupling  spin,  fractional  occupa¬ 
tion  numbers,  and  symmetry  in  a  complete  approximate  LDF  theory22,27,28 
is  unfinished. 

This  work  was  supported  by  the  Office  of  Naval  Research  and  through  a 
grant  of  computer  time  by  the  Research  Advisory  Committee  of  the  Naval 
Research  Laboratory. 
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